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ABSTRACT 

Modam crop growth sioulation and crop 
condition assoaanont nodcls require daily 
input of naxiaua teuperature, oiniMm teu- 
perature, precipitation, and aolac radiation 
data. Gridded spatial eatiBStea of these 
variables arc prepared for agricultural use 
froa World Meteorological Organisation 
surface reports and enhanced by polar 
orbiting satellites. If sufficiently 
accurate, there data which are made available 
fron an operational system soon after they 
are prepared, nay reduce the cost of reliably 
cstiaating world crop conditions. The 
gridded estinates were compared with daily 
netaorological data aeasured at various 
agricultural research facilities across the 
United States to determine their level of 
accuracy. Preliminary results indicate that 
daily maxisnim temperature can be determined 
to within 9.1 degrees Celsius with ninety 
percent confidence. With similar levels of 
confidence, daily minimum temperature can be 
determined to within 6.7 degrees Celsius, 
daily jsolar radiation to within 231.2 
cal/cm min, and daily precipitation to within 
9.7' millimeters. 

1. INTRODUCTION 

This paper outlines a study conducted by 
the Yield Research Branch of the Statistical 
Reporting Service as part of the Agriculture 
and Resources Inventory Surveys Through Aero~ 
space ResK>te Sensing (AgRISTARS) program, a 
joint USDA, NASA, NOAA, USDC, and AID research 
effort to determine the feasibility of inte- 
grating aerospace remote sensing technology 
into existing or future USDA data acquisition 
systems. The purpose of this study was to 
compare the U.S. Air Force (USAF) Agromet data 
to measured daily meteorological data 
collected at various agricultural research 
facilities across the United States. The 
daily data elements evaluated were maximum 
temperature, minimum temperature, total 
precipitation, and solar radiation. Measured 

data are not readily available for evaluation 
of potential evapotranspiration and crop 
growth simulation. If the Agromet data — which 
are routinely available and prepared for 
agricultural use from World Meteorological 
Organisation (WMO) surface reports and 
enhanced by polar orbiting satellites— are 
accurate enough for plant and soil water 
modaling, data collection costs may be 
aignificantly reduced. 


i. BACKGROUND 

Intereat in worldwide crop production and 
economic conditions has grown in recant years. 
In response, the USAF has developed a complex 
model at Global Weather Central (GWC), Offutt 
Air Force Base, Nebraska, to provide daily 
meteorological data specifically tailored for 
agricultural use. Data to provide Agromet 
information comes from the WHO network surface 
reports and polar orbiting satallitaa 
measuring reflectance, radiance, and tempera- 
ture. An automated cljud analysis model 
(30MEPH) estimates the effect of clouds on the 
radiation balance (refs. 1 and 2). 

Through an agreement between the USAF and 
the U. S. DepartsMttt of Agriculture (USDA), 
these Agromet data are available for crop 
condition asiessment and research use on a 
real tine basis. No assessmeac of the accu- 
racy of these lata over a season has been 
done . 

3. DESCRIPTION Of DATA 

The data alements described include only 
those for %«hich an evaluation is being done. 
Agromet uses the 1/8 mesh Af GWC grid on a 
polar stereographic projection, which gives 
approximately 25-nm grid point spacing at 
60 N, and all data elements arc provided on 
this grid point basis. These estimates are 
currently provided for most of the United 
States and many areas of the Northern 
Hemisphere. Preparation of estimates for the 
United States began in June 1981. 

3.1 MAXIMUM AND MINIMUM TEMPERATURES 

Maximum and minimum temperatures are esti- 
mated every 3 hours from satellite temperature 
estimates. The highest and lowest estimates 
are saved for use at the end of the day. 
Surface reports (three hourly temperatures and 
daily reported maximum-minimums ) are used in 
an analysis with the satellite temperature 
estinates to create maximum -minimum tea^era- 
tures for each grid point ir. each geographic 
region. 

3.2 P RECIPITATION 

Precipitatioi reports from surface stations 
are each assigned to the nearest grid point. 
Reported values are accumulated along with 
ettimatad amounts based on weather conditions 
from surface repor ts. Daily accumuiations arc 
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Mdc of tht groottr of tho abovt anounta. 
Chackt ara aade for convactive precipitatic>n, 
and ona-half of tha daily accumulacion it 
apraad to adjacant grid points if no othar 
pracipitation it rtporcad for thota grids. 
Spraading it batad on atciaatad anountt from 
tho 3DKEPH cloud analysis. A ratio of 
raportad to attiaatad precipitation (R/E) is 
calculated for each grid point and spread 
using a linear distance veighting. 

The raportad precipitation is used for those 
•points for which it is available. At othar 
points tht R/E is used to detarmina a value. 
In a few cases the quantitative precipitation 
forecast (QPF) from tha 3DNEPH is used. This 
is generally less than 12 grid points for an 
area. 

3.3 SOLAR RADIATIOX 

Clear sky direct solar radiation is calcu- 
lated from long standing, well-known aqua- 
tions. The clear sky solar radiation is ad- 
justed for cloudiness using 15 cloud layers in 
the 3DNEPH cloud model. A detailed explana- 
tion of Agroset net solar radiation conpu'a- 
tion is given in CSDA EIAC/TN-31-001, March 
1981. 

3.4 OTHER ACROMET DAT A 

Othar data elements produced by Agroaet are 
not being evaluated at this time. A descrip- 
tion of these data can be obtained from 
ETAC/TN-61-001. 

4. DATA COLLECTION 
4.1 AIR FORCE ACaO’lET DATA 

The Agrooet data being used in this study 
ware processed in the Foreign Agricultural 
Service (FAS), Foreign Crop Condition Assess- 
ment Division (FCCAC), Houston, Texas. 

Agromet data are processed continuously 
(every 3 hours) at GVC, Offutt Air Force Base, 
Nebraska. Each day at 2400 CMI the dnily 
Agromer Data Sursary is prepared. Oats from 
these daily reports are assembled every 

Monday into a weekly (Monday throug- Sunday) 
data set, ana a magnetic tape is Air Expressed 
to FCCAD, Houston, Texas. The data are 
processed each Tuesday end are available for 
operational use generally within iO hours 
after preparation at CV>'C. Data for comparison 
to measured date have been extracted from the 
FCCAD disk file for use in statistical 
analysis. 

Grid alesencs for each grid cell in which a 
research location is situated were compared 
against the corresponding measured data 
alamaot for the saiac day. 

4.2 MEAS’.’REO DATA 

The ground-ceasured me : eorolcg i c il data 
vara collected from -routine mcasureser.rs made 
at various agricultural cesearch canters in 
tha United States and were a-'sembled by the 
USDA/Agricultural Resesrch Service (.ARS), Crop 
Syatomi Evaluation Unit at lemple, Texai. 


Daily measuremanta included '■aximum-minimum 
camporaturaa, .irecipitation, and aolar radia- 
tion* 

Both of theae data aata vara furniahad to 
tha Yield Modal Development (YMO) project, 
Houaton, for use in ststistical evaluation of 
tha Agroaat data. 

5. STATISTICAL ANALYSIS 

Analytic and graphical methods ware 
uaad to characterize tha error struetura of 
the Agromet estimators. For each variable tha 
mean, variance, and other atandard acatiacica 
ware computed for differences betwaan the 
Agromet eeasuraments and tha station maaaura- 
menta. The mean and variance of the differ- 
ences provide estimates of the bias and vari- 
ance expected in uaing the Agromet data, 
assuming that the difference estimates the 
error made by using the Agromet data instaad 
of the station data. Tests for month and sta- 
tion effects were examined using a two-factor 
linear model with month and station random 
factors. Similar nonparamctric casts ware 
performed and the results compared with the 
parametric results. Histograms and time ploCt 
vara made and visually examined for each taC 
of differences. Theae graphs provided inaight 
into the expected behavior of Che estimators. 

Since the measured data provides values at 
a specific point and the Agromet data provides 
astimates over a large area (no smaller chan 
25 by 25 no), differences in Che fwo data 
values were expected. Understanding chess 
differences is necessary before using the 
Agromet data for large-area agricultural 
estimation without regard to ground station 
information. 

6. SOME PRELIMINARY RESULTS 

Ona way to characterize the errors made in 
using Che daily Agromet cell etcimates as an 
estimator of the station point measureaents ia 
to fit a general linear model to the daily 
differences between these measurements. Since 
tha station locations are a subset (more or 
less rsodomly located in the major agricul- 
tural land of Che United States) of a much 
larger sec of points for which we wish to make 
inferences the classification variable station 
location was ccr.siiered a r&r.iom effect. Row 
to treat the classification variable month was 
not so obvious. The months themselves were 
not a random sample of all months; the were 
the cu-:lve consecutive months Cor which the 
data were fathered. However, Che effect on 
Che differences was considered s random effect 
becauac we were not iatercated ia comparing 
specific m.)-'hs; we were interested only in 
ascertaining how ouch of the estimation error 
was attributable to monch-eo-month changes in 
the d’Cf»rinc»s. Th- o-.alyjis for Che vari- 
able daily maximum temperature is presanCad as 
an exs--le of the tevi-.nique used and tha 
raiulcs obtained. 

Tne two-f*ctor rando-n-ef fects i:odel was fit 
CO the data. The ccmacnenca of variance 
associated vith the resv';-’!l error, tho conth 
effect, the station effect, and the station 
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«4 aOBth inttractioa wtr« ••tiaatad bjr tkraa 
ataadard aathods (raf 3, pp. 433), and hypo-' 
thataa imra taatad utiod aavaral analegoua 
typa tuas of aquarat to aacartaia if any of 
Cha co ap oa aata of varianca could raaaoaably ha 
coaaldarad to ha aero (rafa 4 and 3)* Tha 
raaulta for taating tha hypothaaia that a 
variaaca coi^oaaBt ia aaro era auaaariaad ia 
Ittla 1; tha raaulta era givaa uaiag typa IV 
aaalo|oua auaa of aquaraa. lha raaulta were 
idaatacal whaa other typa awa of aquara were 
iiaad. Table 2 glvaa the aatiaataa for the 
variaaca coapoaaata for each of tha three 
aathoda aaployad. Tha concluaioa ia that the 
variaaca coapooeat aaaoeiatad with tha claaai- 
ficatioa variable amith ahould be conaidarad 
aero. Following ataadard atatiatical proce- 
dure, tha claaaificacicm variable aoath waa 
droned froa tha analyaia and the reduced 
aodal refitted to the data, tha hypothaaaa 
retaatad, aod tha variance coaponants reea- 
tiaatad. Tba raaulta are presented ia Tablea 
3 and 4. 

lha praliainary concluaion froa the results 
givea ia Tables 1-4 is that tha variance 
eoapoaeats aasociatad with the station loca- 
tioa and tha statim aonth interaction are not 
aaro and chat both eontribuCa significantly to 
the total error aada in using the 4grouet 
daily call Maaianw-tewparature estiwata as an 
•atiaator of the associated atation point 
■aaiauw-taaparature aeasuramenc, Thera ia no 
reason to suapect that the variance co^onent 
aaaoeiatad with the classification variable 
•oath contributes significantly to the overall 
error structure i hence, it is reasonable to 
consider it sero. 

Froa Table 4, tha residual variance, tha 
station variance, and the station*aonth 
variance contribute respectively about 85 
percent, 9 percent, and 7 percent to the total 
variance. This i^liea that about 90 percent 
of Che cine Che daily cell grid aaxinua 
Cenperature estimate is within 9.1 degrees 
Celsius (9.1~1.64 tines the estiaated standard 
deviation) of Che associated point station 
aeaaureaent for daily aaxinua teaperacire ; Che 
coatribuCion of Che residual error to thia 
value is 8.4 degrees Celsuis. 

Siailar results indicate that daily ainiaua 
teaperature can be deterained to within 6.7 
degrees Celsius with 9Q percent confidence. 
With siailar levels of confidence, daily so’ar 
radiation cu be deterained to within 231.2 
(calories/ca ain) and daily precipitation to 
within 9.7 ailliaacera. For daily ainiaua 


teaperature, the residual variaaue, the 
atatioa variance, and the sCaCion*aoath 
variance were estiaated to coatribate 
respectively 7S percent, 12 percent, and 13 
percent to Che total variance. For daily 

solar radiatioa, the residual variance, the 
station variance, and the scatioa^aonch 
variaaca were estiaated to coatributs 
respectively 96 percent, 2 percent, sad 2 
percent to Che total variance. For deity 

precipitation, the residual variaaca and the 
atation variance .were estiaated to contribute 
respectively 99 percent, end 1 pereent to Che 
total variance; the variance uoapo nen t aaoo- 
cieted with the scation*aonCh intaractioa was 
not significant. A eoaplata diaeussion of 
cheae aaalyaea can be found in Ferry 1982 
(ref. 6). 

^Contribution of the Yield Model Develo p a a nt 
(YMD) project within Che Agriculture end 
Sasources Inventory Surveys Through Aaroapece 
loaote Sensing (AgRISTAKS) prograa, a joint 
progren of USOA, USDC, 8ASA, USDl, end AID. 

^AgKISTA?.S, 1050 Bay Area Blvd., HouaCoa, 
Texas 7 /C5B. 

^Grasslands, Soi' and Hater Laboratory, 
teople, Texas 76503 
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TABLE 1.- ANALYSIS OF VARIANCE FOR THE FULL MODEL USING TYPE IV SUMS OF SQUARES TO TEST 
CMq: VAR(NONTH) - 0; Hq: VAR( STATION) - 0; Hg: VAR(STATION*MONTH) • 0] 


SOURCE OF 
VARIATION 

DEGREES OF 
FREEDOM 

SUN OF 
SQUARES 

COMPUTED 

fo 

PROS F > Fq 

EXPECTED MEAN SQUARE 

MONTH 

11 

704.5 

0.67 

m 

o\ * 29.1S!*5^„* ♦ 297.97 

STATION 

10 

9712.5 

10.18 

■ 

c\ * 28.89 * 322.41 oj 

STATION*NONTH 

104 

9921.4 

3.67 

■ 

♦ 29.15 

ERROR 

3551 

92400.9 


■ 
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TABLE 2.< VARIAW:E COMPONENT ESTIMATES FOR THE FULL MODEL 


VARIANCE 

COMPONENT 

VARIANCE COMPONENT ESTIMATION PROCEDURE 

TYPE I SS 

MIVQUEO 

MAXIMUM LIKELIHOTO 

VAR (MONTH) 

- 0.12 

- 0.12 

0.00 

VAR (STATION) 

2.68 

2.69 

2.43 

VAR (STATION*MONTH) 

2.38 

2.39 

2.23 

VAR (ERROR) 

26.02 

26.00 

26.02 


TABLE 3- ANALYSIS OF VARIANCE FOR THE REDUCED MOTEL USING TYPE IV SUMS OF SQUARES TO TEST 
[Hq; VAR(STATION) • 0 ; Hq: VAR(STATI0N*M0NTH) - 0] 


SOURCE OF 
VARIATION 

DEGREES OF 
FREEDOM 

SUM OF 
SQUARES 

COMPUTED 

''o 

PROS F > Fg 

EXPECTED l«AN SOTARES 

STATION 

10 

9897.71 

10.83 


a\ * 28.89 330.39 03 

STATION*MONTH 

104 

10509.79 

3.51 

■ 

9 2 

0 £ * 29«16 

.ERROR 

3551 

92400.94 


WM 

2 

«E 


TABLE 4.. VARIANCE CWPONENT ESTIMATES FOR THE REDUCED MOTEL 


VARIANCE 

COMPONENT 

VARIANCE COMPONENT ESTIMATION PROCEDURE 

TYPE I SS 

MIVQUEO 

MAXIMUM LIKELIHOOD 

ESTIMATE 

PERCENTAGE 

ESTIMATE 

PERCENTAGE 

ESTIMATE 

PERCENTAGE 

VAR (STATION) 

2.71 

8.8 

2.70 

8.7 

2.43 

7.9 

VAR (STATION*MONTH) 

2.24 

7.2 

2.27 

7.3 

2.23 

7.3 

VAR (ERROR) 

26.02 

84.0 

26.00 

84.0 

26.02 

84.8 

TOTAL VARIANCE 

30.97 

100.0 

30.97 1 

100.0 

30.68 

100.0 




















































